Model of Mass Varying Neutrinos in SUSY 
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ABSTRACT 

We discuss the mass varying neutrino scenario in the supersymmetric theory. In the 



> 

^ ] case of the model with the single superfield, one needs the soft SUSY breaking terms or 

t^. 1 the /i term. However, fine-tunings of some parameters are required to be consistent with 
O ' 

■ the cosmological data. In order to avoid the fine-tuning, we discuss the model with two 
O _ 

superfields, which is consistent with the cosmological data. However, it is found that the 

Oh" 

left-handed neutrino mixes with the neutrino of the dark sector maximally. Adding a right- 
handed neutrino, which does not couple to the dark sector, we obtain a favorable model 



in the phenomenology of the neutrino experiments. In this model, the deceleration of the 
cosmological expansion converts to the acceleration near z ~ 0.5. The speed of sound c s 
becomes imaginary if we put wq = —0.9, which corresponds to m° v = 3.17 eV. On the other 
hand, if we take wo = —0.998, which leads to m° = 0.05 eV, c 2 s becomes positive since w 
evolves rapidly near the present epoch in our model. 
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1 Introduction 



One of the most challenging questions in both cosmological and particle physics is the nature 
of the dark energy in the Universe. At the present epoch, the energy density of the Universe 
is dominated by a dark energy component, whose negative pressure causes the expansion of 
the Universe to accelerate. In order to clarify the origin of the dark energy, one has tried to 
understand the connection of the dark energy with particle physics. 

Recently, Fardon, Nelson and Weiner [1] proposed an idea of the mass varying neutrinos 
(MaVaNs), in which the neutrino couples to the dark energy. The variable neutrino mass was 
considered at first in [2], and was discussed for neutrino clouds [3]. However, the renewed 
MaVaNs scenario [1] has tried to make a connection between neutrinos and the dark energy. 
In this scenario, an unknown scalar field which is called "acceleron" is introduced, and then, 
the neutrino mass becomes a dynamical field. 

The acceleron field sits at the instantaneous minimum of its potential, and the cosmic 
expansion only modulates this minimum through changes in the neutrino density. Therefore, 
the neutrino mass is given by the acceleron field and changes with the evolution of the 
Universe. The cosmological parameter w and the dark energy also evolve with the neutrino 
mass. Those evolutions depend on a model of the scalar potential strongly. Typical examples 
of the potential have been discussed by Peccei [4]. 

The MaVaNs scenario leads to interesting phenomenological results. The neutrino oscil- 
lations may be a probe of the dark energy [5, 6]. The baryogenesis [7, 8, 9], the cosmo MSW 
effect of neutrinos [10] and the solar neutrino [11, 12] have been studied in the context of 
this scenario. Cosmological discussions of the scenario are also presented [13, 14, 15, 16, 17]. 

In this paper, we study the MaVaNs scenario in the supersymmetric theory and construct 
models which are consistent with the current cosmological data [18]. We discuss the dark 
energy in the some cases of the superpotential. Then, we present the numerical results for 
the evolution of the neutrino mass and w. In section 2 and 3, we study the superpotential 
with the single superfield and the double superfields, respectively. The section 4 devotes to 
the discussions and summary. 
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2 The single superfield model 

The simplest assumption of the MaVaNs with the supersymmetry is to introduce a single 
chiral superfield A, which is a singlet under the gauge group of the standard model. The 
superfield A couples to the left-handed lepton superfield L. In this framework, we discuss 
three cases of the superpotential. 

2.1 The simplest model W = ^A 3 + m D LA 
We suppose the dark sector with the superpotential 



where A and m D are a coupling constant and a mass parameter, respectively. The scalar 
and spinor components of A are ((p a ,ip a ), and the scalar component <p a is assumed to be 
the acceleron. The second term of the right hand side m D LA in eq.(l) is derived from the 
Yukawa coupling yLAH with y(H) = mp, where H is the Higgs doublet. Assuming the 
vanishing vacuum expectation value of the left-handed slepton, the superpotential of eq.(l) 
is enough to discuss the dark energy, because only the scalar potential of the acceleron and 
the neutrino energy density contribute to the dark energy in the MaVaNs scenario. We omit 
other terms, which does not couple to A, in our superpotential. 
Then, the scalar potential for <p a is given by 



This means that the dark sector interacts with the standard electroweak sector only through 
neutrinos. By solving the eigenvalue equation of the 2x2 mass matrix, 




(1) 



V(<p a ) = \ 2 \<p a t + m 2 D \<p a 



(2) 



We can write down a Lagrangian density from eq.(l), 



C = m D v L i) a + 2A0 a ^ a ^ a . 



(3) 




(4) 



4> a is given in terms of the neutrino mass, 




m 2 D 
2m u 



(5) 
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Using this relation, the scalar potential of eq.(2) is given in terms of the neutrino mass m v 
as follows: 

where, for simplicity, we take the scalar field real. 

In the MaVaNs scenario, there are two constraints on the scalar potential. The first one 
comes from the observation of the Universe, which is that the present dark energy density 
is about 0.7p c , p c being a critical density. Since the dark energy is assumed to be the sum 
of the energy densities of the neutrino and the scalar potential 

Pdark = ^ + y (^( m ")) , ( 7 ) 

the first constraint turns to 

A + V(d>°M))=0.7p e , (8) 

where "0" represents a value at the present epoch and 70% is taken for the dark energy in 
the Universe. 

The second one comes from the fundamental assumption in this scenario, which is that 
/'dark is stationary with respect to variations in the neutrino mass. This assumption is 
represented by 

dm v dm v 

For our purpose it suffices to consider the neutrino mass as a function of the cosmic temper- 
ature [4]. Then the stationary condition eq.(9) turns to [4] 

T3 dF | dV(Umu)) _ Q (1Q) 
<9£ dm u 



where f = m v (T)/T, p v = T 4 F(£) and 

F{ )= i r%W+F. 

7r 2 io ev + 1 K ' 

We have the time evolution of the neutrino mass from the relation of eq.(10). Since the 
stationary condition should be satisfied at the present epoch, the second constraint on the 
scalar potential is 



^8F , dV(Mmu)) 



T 6 ^- + 



<9£ dm u 
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= • (12) 

m v =m°,T=T 



This condition turns to 



= -< , (13) 

m v =m°,T=T 



dm u 

where n° u is the neutrino number density at the present epoch. 

Since neutrinos are supposed to be non-relativistic at the present epoch, p° v = m®n® is 
given and the equation of state becomes 



w° + 1 = — . (14) 

Taking the typical observed value Wq = —0.9, we can fix p° u . Then the neutrino mass is 
obtained by putting the neutrino number density at the present epoch n° = 8.82 x 10~ 13 eV 3 
on pi = mini. Finally, we get m" = 3.17 eV and V((fP a (m? v )) = 2.52 x 10" n eV 4 , where we 
take p^ark = O-^Pc = 2.8 x lCT u eV 4 at the present epoch. The neutrino mass 3.17 eV may 
be large compared with the terrestrial neutrino experimental data. The neutrino mass of 
the 1 eV scale is related with the LSND evidence [19] and will be tested at the MiniBOON 
experiment [20]. On the other hand, putting w = —0.998 we get m" = 0.05 eV, which 
is consistent with the atmospheric neutrino mass scale. Thus, the value of depends on 
wo- In our following analyses, the numerical value of is not so important as far as the 
neutrino is non-relativistic at the present epoch. We take m° = 3.17 eV with w = —0.9 as 
a reference value in the following numerical studies. 

Now, we have two constraints on the potential and its derivative at the present epoch as 
follows: 

V(4>° a (ml)) = 2.52 x 10' 11 eV 4 , (15) 
dV(Mm v )) 



dm u 



= -8.82 x 10" 13 eV 3 . (16) 

,0 



It is found that the gradient of the scalar potential should be negative and very small. These 
constraints on the scalar potential are very severe. By using the potential of eq.(6) in the 
model, we have 

?YpA = i L _!^)7 1 + 4) + ^|L_^)f 1 + 4) . (17) 

om u 4A 2 \ m u J \ mf, J 2A 2 \ m v J \ m z v ) 



Therefore, the scalar potential satisfies the relation 

V{(t> a {m v )) 



dV{(t>a{m v )) 



4 1 i "-D 
m v =mO , T=T 1 ^ (^F 



1 — -P 
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This ratio must be —28.6 from eqs. (15) and (16), however, our input m" = 3.17 eV never 
reproduce this value. One cannot build any models with only one superfield A unless the 
SUSY breaking term or the \i term is added. 

2.2 W = |^4 3 + uidLA with soft breaking terms 

Let us take into account the soft-breaking effect of the supersymmetry. Then the scalar 
potential is given by 

V{<t> a ) = A 2 |0 a | 4 + m 2 D \<p a \ 2 + m 2 |0 a | 2 + Vq , (19) 

where m is the soft-breaking mass and Vq is a constant. The scale of the supersymmetry 
breaking in the standard sector m is supposed to be of order the electroweak scale m — v. 
However, if the dark sector couples to this supersymmetry breaking only via the neutrino, 
radiative corrections give the supersymmetry breaking mass scale of order {m D /v)m. There- 
fore, the soft-breaking mass m in the dark sector is expected to be comparable to mo, which 
is taken to be O(leV). Such a small soft mass of the supersymmetry breaking corresponds 
to the small gravitino mass m 3 / 2 — O(leV), which has been given in the gauge- mediated 
model of ref. [21]. 

The gradient of the potential are given as 

where (f> a is given in terms of m v as in eq.(5). Since we have four free parameters, A, m, Vq 
and rriD, we can adjust parameters to constraints of the potential and its derivative. Putting 
the typical values for two parameters by hand as follows: 

A = 1 , m D = 10 eV , (21) 

with mJJ = 3.17 eV, we have <f>a( m v) — —14.2 eV. Then, m 2 and Vq are fixed by the data of 
eqs. (15) and (16) as follows: 

m 2 = -2 x 14.2 2 - 10 2 + e (eV) 2 , (22) 
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V = 14.2 4 - 14.2 2 e + 0.7p c - p° u (eV) 4 , (23) 

where e = 2\ 2 ((fP a ) 2 + m 2 D +m 2 = -5.67 x 1(T 15 eV 2 . It is remarked that the parameters m 2 
and V are fine-tuned on the order of 10~ 15 eV 2 to guarantee the tiny V(m„) and cH^m^) / <9m^ 
at the present epoch, respectively. Using these parameters, we can get the evolution of the 
neutrino mass and w from the stationary condition and the equation of state, respectively. 
However, since such a case of fine-tunings is not interesting, we do not discuss the case 
furthermore. 

2.3 W = f A 3 + m D LA + § A 2 model 

We consider the model including /x-term as follows: 

W = -A 3 + m D LA + -A 2 , (24) 
3 2 

which leads to the scalar potential as 

V{<p a ) = \\<p 2 a + p<p a \ 2 + m 2 D \<p a \ 2 . (25) 
Taking a Lagrangian density of the form 

£ = m D v L ^ a + (2A0 a + /i)^ a ^a , (26) 

we get (p a in terms of the neutrino mass instead of eq. (5) as follows: 

m v -n m 2 D 

X(pa = ^^-2^- (27) 
The derivative of the scalar potential is given as 

= (l + ^) y t(^a + ^) (2 A0 a + ^) + m 2 ,] . (28) 

It is easily found that there is the parameter set, which satisfies the present data of eqs.(15) 
and (16) as 

[(A0 a + /i)(2A0 a + /i)+m^] ~ l(T 10 eV 2 , m D < 0° ~ /i ~ l(T 3 eV . (29) 

This result indicates the fine-tuning among A0 a , ji and m D on the order of 10~ 7 . If the value 
of m D is much smaller than values of \<j) a and //, numerical solutions are 

a = ±2.24 x i(T 3 eV , n = ^4.48 x 10" 3 eV , \m D \ = 10" 4 eV , (30) 



with A = 1. Such a small value of ~ 10 3 eV may be explained by the suppression of 
M| eV /M planck [22]. 

Since two mass eigenvalues are almost degenerate due to > 2A0 a + /i, the left-handed 
neutrino v L mixes maximally with ip a , which is a kind of sterile neutrinos. Thus, this model 
is unfavored in the phenomenology of neutrino experiments. 

3 The double superfields model 

It is very difficult to build a model with the single superfield without fine-tuning of parameters 
of the model. In this section, we introduce two superfields, A and iV [23], which are singlets 
under the gauge group of the standard model. 

3.1 The simple model W = XANN + m D LA + m' D LN 

It is assumed that the dark sector consists of two chiral superfields A and N, whose scalar 
and spinor components are ((f> a ,ip a ) an d (0 n , respectively, with the superpotential 

W = XANN + m D LA + m D LN , (31) 

where the scalar component 4> a of A is assumed to be the acceleron. The scalar potential is 
given by 

V{<j> a , <t> n ) = A 2 |0„| 4 + 4A 2 |0a0n| 2 + m^ | 2 + m£|0 n | 2 • (32) 

The gradient of this potential is describing as follows: 

= 8A 2 2 a + 2mU a , (33) 

0<pa 

where we assume the scalar component of two chiral superfields to be real. Then, we can 
write a Lagrangian density of the form 

£ = m D v L ^ a + m D v L i) n + \$ a ^ n ^ n + A0„^ a ^ n . (34) 

Therefore the mass matrix in the coupled system of the left-handed neutrino and the dark 
sector is given by 

/ mo rn' D 

m D A0„ | , (35) 
\ m' D \(j) n \(p a 



in the (vl, ip a , ip n ) basis. The eigenvalue equation gives 

_ m 3 v - (\ 2 (j) 2 n + m 2 D + m%)m u - 2\m D m' D (j) n 

<Pa — n o ' V^J 

— m D 

By using this relation, the potential of eq.(32) is given in terms of the neutrino mass m u . 
Putting A = 1 and m' D = 1 eV by hand, other three parameters are fixed by three constraints 
of m°, Vim®) and dV(m u )/dm u \ m ^ =m o as follows: 

0° = -2.42 x l(T 15 eV , |0 n | = 5.02 x 10~ 6 eV , \m D \ = 3.01eV , (37) 

where mass eigenvalues are obtained 

ml = ±3.17eV , -3.01 x 10~ 6 eV . (38) 

We can see that two neutrinos are degenerate in the mass, in other words, vl and u a mix 
maximally. Therefore, this model is also unfabored in the phenomenology of the neutrino 
experiments. 

3.2 Right-handed neutrino 

Towards a realistic model, we introduce a right-handed heavy Majorana neutrino, which is 
assumed to decouple from ip a and ip n . Then, the effective mass matrix in eq.(35) is modified 
in the (ul, ip a , ip n ) basis as follows: 

/ C LL m D m' D \ 

m D A0 n , (39) 
\ m' D X(f) n \<f) a J 

where Cll is the effective mass given by the seesaw mechanism between the left-handed and 
right-handed neutrinos. The eigenvalue equation gives 



= ml - C LL ml - (\ 2 4> 2 n + m 2 D + m%)m v - 2\m D m' D (p n + C LL \ 2 4> 2 n 

ml-m 2 D - C LL m v 



Putting A = 1, m D = O.OleV and m' D = O.leV by hand, other three parameters are fixed 
by the three constraints of m°, V(m°) and dV (m v ) / dm v \ mv=m o as follows: 

0° = -4.38 x 10~ 12 eV , 0„ = ±5.02 x 10" 5 eV , C LL = 3.17 eV , (41) 
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Figure 1: Plot of the scaled neutrino mass versus the redshift z 



where mass eigenvalues are obtained 



rn 



° = 3.17eV, -3.17 x 10" 3 eV, -9.18 x 1(T 6 eV , 



(42) 



where 3.17 eV is the mass of the active neutrino, and other ones are for sterile neutrinos. 
Actually, the mixing between the active neutrino and sterile ones are tiny. 

The evolution of the neutrino mass is given by using the stationary condition of eq.(10). 
We show the scaled neutrino mass m v /m° v versus the redshift z = T/T — 1 in Fig.l, because 
the absolute neutrino mass is not important as far as the neutrino is non-relativistic at 
the present epoch. As seen in Fig.l, the neutrino mass evolves only 0.1%. This weak <fi a 
dependence of the neutrino mass is understandable in the approximate mass formula: 

where the constant term Cll — 3.17 eV dominates the neutrino mass and the <p a dependence 
is suppressed on the order of m|,/C£ L . 

Once the evolution of m v is given, one can calculate the equation of state parameter w 
as follows: 



„ m n + 
m v ~ Cll + — — 



^£(1 + ^*.) 

^LL V ^LL / 



(43) 



w + 1 



4 - MO 



V(m v ) 
T 4 F(0 



(44) 
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Figure 2: Plot of the equation of state parameter w versus z in the region of (a) z — 
100000 and (b) z = ~ 10. 
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Figure 3: Plot of the energy density of the dark energy and of the matter in the unit of p c 
versus z, where the solid line and the dashed line correspond to the dark energy and the 
matter, respectively 



where 



h(0 



(45) 



F(0 d£ 

The evolution of w versus z is shown in Fig. 2. In order to see the behavior of w near the 
present epoch, we also plot w at z = ~ 10. It is noticed that w evolves rapidly near the 
present epoch. 

The evolution of the dark energy in the unit of p c is shown in Fig. 3, in which the evolution 
of the matter is presented in comparison. It is found that the matter dominates the energy 
density of the Universe at z > 1. 
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Figure 4: Acceleration a/ a versus z 

In order to see when the acceleration of the cosmological expansion begun, we calculate 
the acceleration a/ a in the Friedmann equation; 

a = ~~Y~^ PM + ^ W + 1 ' )Pd ark] ' ( 46 ) 
where pu is the matter density and the contribution of radiation is neglected since we 
consider the epoch of z — ~ 1. As seen in Fig. 4, the deceleration of the cosmological 
expansion converts to the acceleration near z = 0.5 in this model. This result is different 
from the one in the power-law or exponential potential discussed by Peccei [4] , in which the 
conversion from the deceleration to the acceleration is predicted near z = 5 ~ 7. 

4 Discussions and Summary 

In our work, we have presented numerical results in the case of the non-relativistic neutrino 
at the present epoch. However, it was remarked that the speed of sound, c s , which is given 
as [15] 

c 2 s =w + j— p dark , (47) 
^dark 

becomes imaginary in the non-relativistic limit at the present and then the Universe cease to 
accelerate. Actually, c 2 s is negative in our potential if we put wo = —0.9, which corresponds 
to m° = 3.17 eV. On the other hand, if we take w = —0.998, which leads to m° = 0.05 eV, 
we get positive c 2 s since the time evolution of Pdark becomes slower in the case of smaller 
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m v near the present epoch. Therefore, the atmospheric mass scale of the neutrino mass 
m° = 0.05 eV may be favored. 

We have not discussed the quantum corrections to the scalar potential. These correc- 
tions were discussed in ref.fl], in which it was remarked that the neutrino mass should be 
lower than O(leV) although these are model dependent. The quantum corrections will be 
investigated carefully in the coupled system of the neutrino and the acceleron [24]. 

We have discussed the MaVaNs scenario in the supersymmetric theory and found a model 
which is consistent with the cosmological data. In the case of the model with the single 
superfield, one needs the soft SUSY breaking terms or the \i term. However, fine-tunings of 
some parameters are required to be consistent with the cosmological data. 

In order to avoid these defects, we have discussed the model with two superfields, which 
is consistent with the cosmological data. However, the left-handed neutrino mixes with the 
neutrino of the dark sector maximally in this case. 

Adding a right-handed neutrino, which does not couple to the dark sector, we obtain the 
model, in which the mixing between the left-handed neutrino and the neutrino of the dark 
sector is tiny. This model is the first example of the MaVaNs with the supersymmetry. In 
our model, the deceleration of the cosmological expansion converts to the acceleration near 
z = 0.5. The related phenomena of our scenario and the extension to the three families of 
the active neutrinos will be discussed elsewhere. 

After this paper was submitted to the journal, the paper in ref. [25] appeared, which 
presents a different supersymmetric model of MaVaNs. 
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